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Abstract 

Background: Meclianical transmission of tlie major livestock pathogen Trypanosoma vivax by other biting flies than 
tsetse allows its spread from Africa to the New World. Genetic studies are restricted to a small number of isolates 
and based on molecular markers that evolve too slowly to resolve the relationships between American and West 
African populations and, thus, unable us to uncover the recent history of T. vivax in the New World. 

Methods: T. vivax genetic diversity, population structure and the source of outbreaks was investigated through the 
microsatellite multiloci (7 loci) genotype (MLGs) analysis in South America (47isolates from Brazil, Venezuela and 
French Guiana) and West Africa (12 isolates from The Gambia, Burkina Faso, Ghana, Benin and Nigeria). 
Relationships among MLGs were explored using phylogenetic, principal component and STRUCTURE analyses. 

Results: Although closely phylogenetically related, for the first time, genetic differences were detected between 
T. vivax isolates from South America (1 1 genotypes/47 isolates) and West Africa (12 genotypes/12 isolates) with no 
MLGs in common. Diversity was far greater across West Africa than in South America, where genotypes from Brazil 
(MLGl-6), Venezuela (MLG7-10) and French Guiana (MLGll) shared similar but not identical allele composition. No 
MLG was exclusive to asymptomatic (endemic areas) or sick (outbreaks in non-endemic areas) animals, but only 
MLGsl, 2 and 3 were responsible for severe haematological and neurological disorders. 

Conclusions: Our results revealed closely related genotypes of T. vivax in Brazil and Venezuela, regardless of 
endemicity and clinical conditions of the infected livestock. The MLGs analysis from T. vivax across SA and WA 
support clonal propagation, and is consistent with the hypothesis that the SA populations examined here derived 
from common ancestors recently introduced from West Africa. The molecular markers defined here are valuable to 
assess the genetic diversity, to track the source and dispersion of outbreaks, and to explore the epidemiological 
and pathological significance of T. vivax genotypes. 
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Background 

Animal trypanosomosis (Nagana) caused by Trypanosoma 
vivax can be a highly debilitating disease in African and 
South American livestock [1,2]. In Africa, T, vivax is 
highly prevalent in both tsetse-infested and tsetse-free re- 
gions. It is considered an important pathogen in Burkina 
Faso [3,4], Ghana, Zambia [5], Nigeria [6], Uganda [7], 
Ethiopia [8,9], Sudan [10] and Cameroon [11]. The cyc- 
lical transmission of T, vivax is limited to tsetse flies; 
mechanical transmission by other biting flies allows 
T, vivax to spread in some tsetse-free African regions 
and to Central and South America, where it is dissemi- 
nated by tabanids and stomoxes [1,2]. In South America, 
the transplacental transmission of T, vivax also plays an 
important role in its epidemiology [12,13]. 

In Africa, bovids and suids are hosts of T, vivax, and 
this species can be pathogenic in equines, camels, cattle, 
goats and sheep, whereas wild ungulates serve as reser- 
voirs [14-19]. In South America, wild reservoirs are un- 
known and T. vivax can be pathogenic to cattle, sheep, 
goats and horses [12,20-23]. The main manifestation of 
acute T. vivax disease is devastating anaemia and com- 
promised cardiac function [14,15,22,24,25], followed by 
the invasion of the central nervous system [20,23,26,27], 
abortion, stillbirth, and testicular and ovarian damage 
[12,13,28]. Disease severity and particular clinical signs 
have been associated with geography, prior infections, 
health conditions, and livestock species and breeds. In 
general. West African (WA) isolates are more pathogenic 
to livestock than East African (EA) T. vivax isolates, but 
wasting disease with haemorrhagic syndromes has been 
reported in Kenya and Uganda [14,15,24]. Outbreaks of 
acute haematological and neurological disorders with 
high mortality have been reported affecting cattle, 
goats, sheep and horses throughout non-endemic Brazil- 
ian regions [20-22,29] . In South American (SA) regions of 
enzootic stability (Amazonian lowlands, Venezuelan 
Llanos and Brazilian wetland of the Pantanal), infections 
are mostly asymptomatic in cattle, buffaloes and sheep, all 
showing low parasitaemias [1,30-33]. 

T. vivax is endemic in many countries in Central and 
South America. The first reports of T. vivax in the New 
World were in French Guiana (1919), Venezuela (1920), 
the Caribbean Guadalupe and Martinique (1926 and 
1929), and Colombia (1931), which are former French 
and Spanish colonies. In Brazil, a former Portuguese col- 
ony, T. vivax was first recorded in cattle (1946) and buffa- 
loes (1972) in Amazonia [1,2]. Based on historical livestock 
introduction and limited parasite genetic evidence, it has 
been suggested that T. vivax was introduced into the 
Americas via West African cattle [1,2,30,34-37]. Cattle, 
horses, sheep, goats, donkeys and pigs were first brought 
to the Americas on the voyage of Columbus in 1493, and 
for centuries the transport of Iberian and African livestock 



to the Americas followed the routes of the African slave 
trade [38-40]. Therefore, cows, goats, sheep and equines 
brought by the colonisers could be responsible for the 
introduction of T. vivax into the Americas at different 
times and places. 

Early studies comparing African and American T. vivax 
isolates using molecular markers showed a relevant simi- 
larity between SA (Colombian) and WA isolates [36,37]. 
The close relationship between SA and WA T. vivax was 
corroborated by phylogenetic analyses of cattle isolates 
from Brazil (3 isolates), Venezuela (one isolate). West Af- 
rica (Y486 from Nigeria) and East Africa (IL3905 from 
Kenya) using Spliced Leader [30], SSU and ITS rDNA 
[34], and Cathepsin L-like [35] sequences. Previous studies 
evidenced high divergence separating SA/WA from EA 
parasites, and also revealed substantial divergence among 
EA isolates from Kenya and Mozambique [30,34,35]. 
Highly divergent isolates from Tanzania (EA) were re- 
ported from tsetse flies by comparing gGAPDH sequences 
[41,42], and also in wild animals through ITS rDNA ana- 
lysis [19]. Cattle isolates from Ghana clustered with SA/ 
WA isolates, while Zambian isolates were more related to 
Kenyan T. vivax [5]. Genetic studies have revealed more 
relevant genetic diversity in EA T. vivax compared to pop- 
ulations in WA [17,19,30,34,35,41-43]. Unfortunately, the 
use of different molecular markers prevented a global 
comparison of data from all previous studies. 

The genetic studies on T. vivax were mostly based on 
molecular markers that evolve too slowly to resolve the re- 
lationships between SA and WA populations and, hence, 
were unable to uncover the recent history of this parasite 
in the New World. Microsatellite multiloci genotype 
(MLG) analysis can reveal cryptic genetic diversity, popu- 
lation structure and the origin of parasites, as have been 
shown for T. brucei spp. [44-47] and T, congolense [48,49] . 
The most comprehensive genetic study of T, vivax by 
MLG analysis was restricted to isolates from donkeys in 
The Gambia and results suggested a clonal population 
[50]. More considerable MLG polymorphisms were dem- 
onstrated in Cameroon [51] and Uganda [52], despite the 
few isolates examined. 

In this study, we analysed polymorphisms in 7 MST 
loci in isolates from across South America (39 from Brazil, 
7 from Venezuela and one from French Guiana) and West 
Africa (12 isolates from The Gambia, Burkina Faso, 
Ghana, Benin and Nigeria) aiming to assess the genetic 
repertoire and phylogenetic relationships at continental 
and intercontinental levels and, hence, to understand the 
introduction and dispersion of T, vivax in South America. 
The comparison of isolates from asymptomatic livestock 
living in areas of enzootic stability and isolates from sick 
animals from outbreaks, exhibiting a range of haemato- 
logical and neurological signs and several fatal cases, as 
well as repeat sampling from the same areas, allowed us 
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to examine potential links between genotype and disease, 
outbreaks, host species and virulence in the context of 
spatial-temporal changes. 

Methods 

T. vivax samples from endemic settings and outbreaks 

The identification of T. vivax in blood samples was per- 
formed using a T, vivax-specific PGR assay (TviCATL-PCR) 
based on Cathepsin-L gene [35]. The South American iso- 
lates were from blood samples collected at widely distrib- 
uted locations in Brazil and Venezuela (Figure 1) from 
cows, buffaloes, sheep and horses. In addition, one sheep 
isolate from French Guiana was included (Table 1). All SA 
isolates, including their host species and geographic ori- 
gin, and clinical signs of infected livestock are detailed in 
an additional table [see Additional file 1]. 

A total of 47 South American isolates were genotyped 
by MST analyses: 39 from Brazil, 7 from Venezuela, and 
one from French Guiana. Furthermore, 12 samples from 
West Africa and two from East Africa were also analysed. 
T. vivax isolates from animals infected during different 
Brazilian outbreaks (Figure 1) were compared. Isolates 
from very sick sheep were from an outbreak in the north- 
east (state of Paraiba) that affected a flock and caused 



severe haematological and neurological alterations 
along with several deaths [20]; other outbreaks in cows 
and sheep occurred previously in the same region [29]. 
Two isolates were from a cow and a horse that showed 
haematological and neurological disturbances in sequen- 
tial outbreaks in the southern region (state of Rio Grande 
do Sul) [21]. Two isolates were from an outbreak in the 
south-eastern region (Sao Paulo) from cattle that exhib- 
ited severe haematological changes [22]. The Venezuelan 
isolates were from cattle, buffalo and sheep from the 
endemic Llanos region (Figure 1), where the infected 
animals are generally asymptomatic but can sporadic- 
ally show moderate parasitaemia and anaemia [32,33]. 
The isolates from West Africa included in this study were 
all from low parasitemic cattle showing mild anemia and 
were collected in T, vivax endemic settings in Burkina 
Faso, The Gambia, Ghana, Benin and Nigeria. Isolates 
from Mozambique were obtained from cattle and a nyala 
antelope [17]. The field-collected samples were designated 
as primary samples, and the laboratory samples represent 
isolates expanded in experimental animals (Table 1). 
T, vivax from Africa, including host and geographic ori- 
gin, and clinical signs of infected animals are detailed in 
an additional table [see Additional file 2]. 




Figure 1 Geographical origin of Trypanosoma vivax isolates from South America and Africa and microsatellite genotype (MLGs) 
distribution. Map showing the different localities of T. vivox isolates used in this study, including samples from endemic regions and from 
outbreaks (*) in non-endemic areas. Brazilian States: PA, Para; PB, Paraiba; MS, Mato Grosso do Sul; SP, Sao Paulo; RS, Rio Grande do Sul. Coloured 
symbols (boxes) correspond to individual MLG or clusters of MLGs defined with the analysis of 7 microsatellite loci (Figures 2 and 3). 

V J 
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Table 1 Geographical and host origin, clinical manifestation and microsatellite loci genotyping data from Trypanosoma 
vivax isolates included in this study 



Geographical origin 



Host species 
(N° of isolates) 



Clinical 
manifestation 



MLG 



Genotypic diversity 
(N° of MLGs/isolates) 



Allelic composition 
(N° of different alleles) 



South America (SA) 
Brazil 

Endemic area 
Non-endemic area 



(outbreal<) 
(outbreal<) 
(outbreal<) 
Total 

Venezuela 

Endemic area 



Total 

French Guiana 

Endemic area 
Total SA 

West Africa (WA) 

Burl<ina Paso 

Gliana 

Benin 

Nigeria 

Tlie Gambia 

total WA 

East Africa (EA) 

Mozambique 

Total EA 

Total SA + WA + EA 



Buffalo (17) 
Cow (4) 
Buffalo (1) 
Cow (1) 
Sheep (2) 
Cow (2) 
Sheep (11) 
Horse (1) 
(39) 

Cow (3) 
Buffalo (2) 
Sheep (2) 
(7) 

Sheep (1) 
(47) 

Cow (4) 
Cow (4) 
Cow (2) 
Cow (1) 
Cow (1) 
(12) 

Cow (1) 
Nyala (1) 
(2) 
61 



Asymptomatic^ 

Asymptomatic^ 

asymptomatic^ 

Asymptomatic^ 

Asymptomatic'^ 

Severe/lethal^ 

Severe/lethal^ 



Asymptomatic^ 
Asymptomatic^ 
Asymptomatic^ 



Symptomatic^ 



Symptomatic^ 
Symptomatic^ 
Symptomatic^ 
Symptomatic^ 
Symptomatic^ 



Symptomatic^ 
ND^ 



1,3,6 
1, 2 
5 
5 

4,5 
1 

2 
3 



7,9,10 
8,10 



12,13,14,15 

16,17,18,19 

20,21 

22 

23 



24 
25 



0.15 (6/39) 



0.57 (4/7) 

(1/1) 
0.23 (11/47) 

(4/4) 
(4/4) 
(2/2) 
(1/1) 
(1/1) 
1.0 (12/12) 

(1/1) 
(1/1) 
1 .0 (2/2) 
0.41 (25/61) 



16 



15 

12 
21 



41 



17 



a, normal PCV and very low 
parasitemia. Details of each 



parasitemia; b, high parasitemia severe anemia (low PCV values) and nervous signs; c, low PCV values and low parasitemia; d, high 
isolate are in Additional files 1 and 2. 



Ethical approval 

The handling of livestock was performed in strict ac- 
cordance with good animal practice as defined by the 
World Organization for Animal Health guidelines and 
approved by Veterinary Scientific Boards of the Centre 
International de Recherche Developpement sur TElevage en 
zone Sub-humide in Burkina Faso, Universidade Eduardo 
Mondlane in Mozambique, and Brazilian Universities that 
participated in this study. The whole project was con- 
ducted in strict accordance with the recommendations of 
the Brazilian National Council of Animal Experimentation 



(http://wwwxobea.org.br/) and approved by the Animal 
Experimentation Ethics Committee from the Institute of 
Biomedical Center, University of Sao Paulo, Brazil (CEP- 
ICB n° 317/09). 

Microsatellite markers and analyses 

Data from the T. vivax Y486 Genome at the Sanger 
Institute (http://www.sanger.ac.uk/Projects/T_vivax/) were 
used for searching MST through the Microsatellite Re- 
peats Finder (http://tandem.bu.edu/trf/trf.html) program. 
MST loci were selected from 5 different scaffolds reducing 
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physical linkages. We designed primer pairs for 14 loci, 
and one primer from each locus was labelled with FAM. 
An additional table shows the sequences of all primers, 
genes and genome locations, and MST motifs of all loci 
examined [see Additional file 3]. 

The 14 primer pairs were initially tested using purified 
DNA from the reference Y486 T. vivax, two further iso- 
lates from Brazil and Mozambique, and other species 
(r. b, bruceU T, evansi and T, congolense). Seven of the 
14 primer pairs (MST loci 4, 7, 8, 10, 11, 13 and 15) 
were specific for T. vivax. The PGR amplifications were 
performed in a 25 \A reaction mixture consisting of -20 ng 
of DNA, 100 pmol of each primer, 200 mM of each dNTP, 
10 mM Tris-HCl (pH 8.3), 3.0 mM MgCL2, 7.5% (v/v) di- 
methyl sulphoxide, 0.1 mg/ml bovine serum albumin and 
1.0 U Taq DNA polymerase. The amplification conditions 
were as follows: initial denaturation at 95°C for 3 min 
followed by 30 cycles of 30 s at 95°C, 30 s at the specific 
annealing temperature for each marker, 1 min at 72°C, 
and a final extension at 72°C for 5 min. The annealing 
temperatures were: 55°C (MST 4, 8 and 11); 58°C (MST 
13 and 15); or 60°C (MST 7 and 10). The allele sizes for 
each locus were determined using a capillary-based se- 
quencer and the Gene Mapper® software with Gene Scan 
500-ROX size standards (Applied Biosystems). The indi- 
vidual peaks defined each allele, and the data set from the 
7 loci defined each MLG. 

Allele frequencies and estimates of genetic variation 
within populations (average numbers of alleles per locus, 
allelic richness and the means of the expected, total ex- 
pected and observed heterozygosity) were calculated 
using ARLEQUIN 3.5 [53]. Genotypic diversity was esti- 
mated as the number of different MLGs divided by the total 
number of isolates. Conformation to Hardy- Weinberg 
equilibrium, as a test of the non-random associations of 
alleles within diploid individuals, and the linkage disequi- 
librium between all pairs of loci, as a test of the non- 
random association of alleles at different loci, were also 
determined in ARLEQUIN 3.5 [53]. The Fixation Index 
(Fit), as a measure of an overall inbreeding coefficient, 
was determined in GenAlEx 6 [54]. Fit values range 
from -1 to 1, where values close to zero are expected 
under random mating, substantial positive values indi- 
cate inbreeding, and negative values indicate an ex- 
cess of heterozygosity. 

To examine the relationships between SA and WA 
T, vivax, the pairwise measure of shared allele distances 
from the microsatellite dataset was calculated using the 
program POPULATIONS vl.2.30 beta [55], and dendro- 
grams based on MLGs were constructed using the DAS, 
shared allele distance [56] and the neighbour-joining (NJ) 
method (bootstrap based on 100 replicates). Principal 
component analyses (PCA) of the MLGs were performed 
in GenAlEx 6 [54]. A Bayesian clustering approach as 



implemented in the program STRUCTURE v2.3.3 [57] 
was employed to estimate the number of genetically differ- 
entiated clusters (/<) within the data set. Log-likelihood 
values for each value of K (ranging from 1 to 16) were 
evaluated from all MLGs by running the STRUCTURE 
program with 300,000 repetitions for three replicates 
(burn-in = 100,000 iterations), and the most likely value of 
K was assessed by the method of Evanno et al (2005) 
[58]. STRUCTURE analysis was performed as described 
previously for T, brucei ssp. MLG studies [45-47]. 

Results 

Microsatellite multilocus genotyping and relationships of 
r. vivax in South America and West Africa 

The MLG analysis of T, vivax populations included 47 
SA and 12 WA isolates, most genotyped directly from 
blood samples, thereby avoiding parasite selection by in- 
oculation in experimental animals. Laboratory isolates 
submitted to successive passages in animals were in- 
cluded for comparison (Table 1, Additional file 1). The 
total number of MLGs found was 25, consisting of 11 
MLGs in 47 SA samples (MLGs 1-11) and 12 MLGs in 
12 WA samples (MLGs 12-23). Therefore, WA T. vivax 
displayed considerably greater genetic diversity than 
the SA parasites. Among the SA isolates, MSTll and 
MST 15 were invariant, the MLGs generally differed at 
only one MST locus, and many samples shared the 
same MLG (Figure 2). An additional file shows results 
for all SA isolates [see Additional file 1]. The WA MLGs 
were defined by unique alleles found in Burkina Paso 
(MLGs 12-15), Ghana (16-19), Benin (20 and 21), Nigeria 
(22) and The Gambia (23). Although no MLG was shared 
between SA and WA T, vivax, most alleles were shared by 
the two populations. The SA isolates showed high homo- 
geneity, with small variability restricted to four loci. While 
12 alleles were shared between WA and SA, only 4 alleles 
were shared between EA and SA T, vivax. The two EA 
isolates included in the analyses were assigned to two 
unique MLGs (24 and 25) clearly separated from both SA 
and WA MLGs (Table 1; Figure 2). 

To illustrate the MLG repertoire and allelic compos- 
ition from SA isolates, 16 samples representing all allele 
profiles (MLGs) identified across the host species and 
geographical range were compared with the African iso- 
lates (Figure 2). An additional file shows the data from 
the remaining SA isolates [see Additional file 1]. Eleven 
MLGs (1-11) were identified in SA and the closely re- 
lated MLGs 1-6 were detected exclusively in Brazil. 
MLGs 1-2 were the most prevalent (-79.5%), probably 
due to many samples from Amazonia and Semiarid. 
MLG 1 was found in 70% of the Amazonian asymptomatic 
buffaloes, which also harboured MLGs 3 and 6, whereas 
MLG 2 was detected in asymptomatic cattle in this region. 
No MLG was exclusive of sick animals. MLG 1 was found 
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T. vivax 



MST loci and allele sizes 





isolates 




MST4 


MST7 


MST8 


MST10 


MST11 


IV1ST13 


MST15 




TviBrPA20 


E 


236 


236 


226 


242 


178 


184 


115 


115 


147 


155 


144 


146 


172 


18^ 




TviBrSPI 


NE 


236 


236 


226 


242 


178 


184 


115 


115 


147 


155 


144 


146 


172 


186 




TviBrPB27 


NE 


236 


236 


226 


242 


178 


184 


115 


115 


147 


155 


144 


148 


172 


186 




TvlBrPASO 


E 


236 


236 


226 


242 


178 


184 


115 


115 


147 


155 


144 


148 


172 


186 


Brazil 
SA 


TviBrPA48 
TviBrRSI 


E 
NE 


234 
234 


234 
234 


228 
228 


242 
242 


178 
178 


184 
184 


113 
113 


113 
113 


147 
147 


155 
155 


144 
144 


146 
146 


172 
172 


186 
186 




TvlBrPBSO 


NE 


236 


236 


228 


246 


178 


184 


115 


115 


147 


155 


144 


148 


172 


186 




TviBrPBSI 


NE 


234 


234 


228 


242 


178 


184 


115 


115 


147 


155 


144 


148 


172 


186 




TviBrPA46 


E 


234 


234 


228 


242 


178 


184 


115 


115 


147 


155 


144 


146 


172 


186 




TviVzApl 


E 


234 


234 


242 


242 


180 


184 


115 


115 


147 


155 


144 


146 


172 


186 




TviVzAp2 


E 


236 


236 


228 


242 


180 


184 


115 


115 


147 


155 


144 


146 


172 


186 


Venezuela 
SA 


TviVzGul 
TviVzAnzl 


E 
E 


236 
234 


236 
234 


228 
228 


242 
242 


180 
174 


184 
184 


115 
115 


115 
117 


147 
147 


155 
155 


144 
144 


146 
146 


172 
172 


186 
186 




TviVzAnz3 


E 


234 


234 


228 


242 


180 


184 


115 


115 


147 


155 


144 


146 


172 


186 


French 
Guiana 


TviVzCoj 
_ TviGuyane 


E 
E 


234 
236 


234 
236 


228 
226 


242 
238 


180 
178 


184 
184 


115 
115 


115 
115 


147 
147 


155 
155 


144 
144 


146 
146 


172 
172 


186 

isJ 


~ TviBfL445 


E 


236 


246 


228 


256 


178 


182 


115 


115 


149 


149 


142 


156 


184 




Burkina 
Faso 
WA 


TviBfMatorkouE 


236 


246 


228 


260 


182 


182 


111 


115 


187 


187 . 






^84 




TviBfFolonzo E 


236 


236 


228 


264 


182 


182 


109 


109 


151 


i 


1 


1 








TvlBfMene 


E 


236 


236 


228 


258 


182 


182 


109 


111 


149 


187 


146 


154 


184 






~ TviKCA19J56 E 


234 


234 


228 


262 


184 


184 


115 


115 


149 


187 




152 


184 




Ghana 


TviDere091 


E 


236 


236 


226 


266 


178 


182 


107 


115 


151 


i 




i 






WA 


TviKC92J28 


E 


234 


246 


226 


260 


182 


184 


m 


109 


149 








184 






_ TviKang92 


E 


232 


232 


228 


228 


182 


184 


107 


135 


149 


175 


156 


158 


184 




Benin 


~ TviBanI 


E 


236 


246 


226 


256 


182 


182 


109 


117 


149 


18^ 


144 








WA 


_ TviBanI. 2 


E 


236 


236 


226 


256 


182 


182 


109 


117 


149 


185 


144 


146 


184 


184 


Nigeria 
WA 

The Gambia 
WA 


I IL700 
Gambia 


E 
E 


236 
234 


236 
234 


228 
228 


256 
248 


182 
182 


184 
184 


109 
109 


127 
111 


163 
163 


183 
183 


142 
142 


144 
144 


184 
184 


186 

,s| 


Mozambique 


" TviMzCb12 


E 


236 


236 


226 


240 


166 


166 


109 


109 


153 


153 


144 


144 


164 


164 


EA 


_ TvlMzNy 


E 


232 


246 


228 


228 


154 


154 


105 


123 


135 


135 


124 


124 


160 


160 



-MLG 



9 

10 

11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 

Geographyc 
distribution 

Figure 2 Microsatellite genotypes (MLGs) of Trypanosoma vivax. The T. vivox isolates included in this figure are representatives of all MLGs 
found in isolates from South American (SA) endemic (E) and non-endemic (NE) areas (16 SA isolates); plus 12 isolates from West Africa (WA) and 
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in fatally infected cows from an outbreak in the south- 
eastern region [22], MLG 2 was detected in a sheep out- 
break of fatal infection with nervous signs in north-eastern 
Brazil [20] where, one year later, MLG 5 was detected in 
asymptomatic sheep, cattle and buffalo, and MLG 4 in 
asymptomatic sheep. MLG 3 was found in a horse with 
fatal infection in the southern region where cattle were re- 
ported infected with T. vivax three years before [21]. 

To represent the relationships between T, vivax isolates, 
we constructed a NJ dendrogram based on pairwise 
distances between the MLGs (Figure 3). This analysis 
showed the separation between the WA and SA isolates. 
In addition, most Brazilian isolates sharing highly similar 
MLGs clustered together, whereas the Venezuelan isolates 
clustered closely with some of the Brazilian isolates. 



forming a heterogeneous cluster of isolates assigned to dif- 
ferent MLGs. The WA isolates formed a highly heteroge- 
neous cluster; only two samples from Benin and three out 
of the four samples from Burkina Faso grouped together. 
The dendrogram was consistent with the clustering 
analyses (Figures 4 and 5), even though most nodes 
had low bootstrap support, most likely due to the low 
number of isolates and loci examined. 

We further assessed the relative genetic differentiation 
between the clusters using PC A (Figure 4), which also 
allowed for the visualisation of the three main clusters 
evidenced by the inferred dendrogram. PGA revealed a 
cluster of isolates exclusively from Brazil, with other 
clusters formed by isolates from Venezuela and some 
isolates from Brazil; clustering in different quadrants 
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Figure 3 Phylogenetic relationships of Trypanosoma vivax from South America and Africa based on microsatellite genotypes (MLGs). 

Neighbour-joining dendrogram based on sliared allele distance; bootstrap values are shown for the major nodes. Symbols represent the isolates 
and respective MLGs. Colours represent the following groups: a) formed exclusively by homogeneous Brazilian isolates (red) and one isolate from 
French Guiana (red with yellow borders), b) moderately heterogeneous Venezuelan (yellow) and Brazilian (yellow/red) isolates, c) highly heterogeneous 
isolates from West (green) and East (blue/green) Africa. Codes of the isolates in South American (TviBr = Brazil, TviVz = Venezuela, 
TviGuyane = French Guiana), West African (TviBf= Burkina Faso, TviK^viDere = Ghana, TviBan = Benin, IL700 = Nigeria, Gambia = The Gambia) 
and East African (JviMz = Mozambique). 



indicating sub-structuring within SA populations. The 
last cluster consisted exclusively of WA isolates distributed 
in two quadrants (Figure 4). Using the method of Evanno 
et al [58], the inferred most-likely number of genetically 
distinct clusters (K) in the STRUCTURE analysis was 
K = 2, corresponding to very well separated SA and WA 
T. vivax populations. However, additional groups were 



present in K > 3, with the split of all the Venezuelan and 
some Brazilian isolates (Figure 5), in agreement with both 
the NJ dendrogram and PCA analyses (Figures 3 and 4). 

Population structure of 7. vivax in South America 

The population structure of T, vivax in South America 
was examined by analysing genotypic diversity among 
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Figure 4 Principal component analysis (PCA) of Trypanosoma vivax populations. PCA analysis of 61 isolates from South Annerica (47), West 
Africa (12) and East Africa (2) based on data from 7 microsatellite multiloci genotyping (MLG). Multidimensional scaling plots of genetic distances: 
Principal coordinate 1 separates South American and African populations besides evidencing the genetic diversity within both populations, 
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TviBan = Benin, IL700 = Nigeria, Gambia = The Gambia) and East Africa (TviMz = Mozambique). 



isolates, multilocus standardised index of association, 
and inter-population differentiation. All 7 MST markers 
were polymorphic in South America; because none of 
the samples displayed more than two alleles per single 
locus, all samples are diploid and have a single genotype. 
Loci MST-7, MST-8, MST-10 and MST-13 displayed be- 
tween 3 and 5 alleles per locus, whereas MST-4, MST- 
11 and MST- 15 showed two alleles (Figure 2). The allelic 
composition of the SA isolates was 21 for 47 isolates 
examined with 9 unique alleles. Neis unbiased genetic 
diversity (Hg) ranged from 0.1 to 0.6 (average of 0.5) 
assuming neutrality. The results indicate a low genetic 
diversity of T. vivax in South America compared to 
West (41 alleles for 12 isolates with 9 unique alleles) 
and East Africa (17 alleles for only two isolates with 
10 unique alleles) (Table 1). 

The analysis of non-random association of alleles 
revealed a significant deviation (P < 0.05) from Hardy- 
Weinberg equilibrium predictions at all loci, and this 
observation was linked to the global heterozygote excess, 
with most loci at, or close to, heterozygote fixation. As 
expected for heterozygote excess, the Fix values for these 
loci ranged from -0.57 to -1.0. Two loci (MST-4 and 
MST- 10) showed heterozygote deficits, with Fit values of 
1.0 and 0.79, respectively. However, the mean observed 



heterozygosity for SA isolates was higher than expected 
(0.71 vs. 0.47), and the mean Fit value reflects this result 
(-0.31). A significant linkage disequilibrium (P < 0.05) was 
observed in the majority of loci combinations. 

The results support a clonal population structure of 
T, vivax in South America. The limited number of isolates 
prevented statistical analysis of the African populations. 
Our attempts to include several EA samples in the MLG 
analysis were unsuccessful, most likely due to the large 
genetic divergence within EA genotypes precluding the 
use of primers based on the WA T, vivax Y486 genome. 

Cryptic diversity within South American 7. vivax infecting 
livestock in endemic areas and outbreaks 

The MLG analyses of SA T, vivax carried out in this 
work did not support association with host species or 
clinical manifestation: buffaloes harbour MLGs 1, 3, 5 
and 6; cattle harbour MLGs 1, 2 and 5; and sheep 
harbour MLGs 2, 4 and 5 (Table 1). An additional file 
shows the MLGs designed for all SA isolates character- 
ized in this study [see Additional file 1]. 

MLGs 7-10 were exclusive to Venezuelan samples 
from animals in endemic regions: MLG 8 and MLG 10 
were found in buffaloes, MLG 7, 9 and 10 in cattle, and 
MLG 8 in sheep. The genotype MLG 8 was shared 
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Figure 5 Genetic analysis of Trypanosoma vivax isolates from South America and Africa. Bayesian clustering (STRUCTURE analysis plots) of 
61 isolates based on 7 microsatellite loci are shown for K = 2 (two main clusters corresponding to South American and African populations) and 
K = 4 (partition of South American samples in two main clusters, one formed exclusively by Brazilian isolates and the other comprising all 
Venezuelan and some Brazilian isolates. Each cluster of isolates is represented by a different colour according to geographic region; thin black 
horizontal lines separate the populations. Codes of isolates in South America (TviBr= Brazil, TviVz = Venezuela, TviGuyane = French Guiana), West 
Africa (TviBf= Burkina Paso, TviK = Ghana, TviBan = Benin, IL700 = Nigeria, Gambia = The Gambia) and East Africa (TviMz = Mozambique). 



between buffalo and sheep from neighbouring regions 
(Apure and Guarico), where animal interchange is in- 
tense, whereas more divergent genotypes were found in 
cattle from a more distant region (Anzoategui) (Figure 2). 
The results suggested that some spatial sub-structuring 
separated Brazilian and Venezuelan isolates. Neverthe- 
less, some samples from Brazil were more closely related 
to those from Venezuela than to other Brazilian samples 
in the NJ dendrogram, PCA and STRUCTURE analyses. 
The only isolate from French Guiana included in this 
work was assigned to an exclusive genotype (MLG 11) 
(Table 1; Figures 2 and 3). 

Discussion 

This study of the population genetic structure of T. vivax 
from South America and West Africa through MLG 
analysis corroborated the close phylogenetic relation- 
ships among SA and WA genotypes previously suggested 
based on a much smaller sample size and limited geo- 
graphical sampling based on different markers such as 
kDNA [37], Spliced Leader [30], SSU and ITS rDNA 



[19,34], Cathepsin L-like [5,35] and gGAPDH [41,42] se- 
quences. The results from MLG analysis of SA and WA 
isolates characterized in this work provides additional 
support to previous hypotheses that have emerged over 
the last decades that T, vivax was introduced from West 
Africa into South America [30,34-37,42]. However, no 
MLG was shared between SA and WA T, vivax, revealing 
that despite closely related SA and WA T, vivax are genet- 
ically distinct. This finding may indicate that mutation 
events occurred in either WA or SA, reflecting discon- 
nected populations, or it may simply reflect the absence in 
this study of the African genotypes that represent the 
source of T, vivax in South America. In addition, the 
MLG analysis revealed small but significant genetic differ- 
ences within SA populations with evidence of some sub- 
structuring: most Brazilian isolates grouped in a single 
cluster, with a second cluster containing all Venezuelan 
isolates along with a few Brazilian isolates from differ- 
ent regions. 

Additionally, our findings suggest, for the first time, that 
genotypes from common ancestry give rise to closely 
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related but genetically different and widespread popula- 
tions in South America. However, despite our analysis 
including samples from former Portuguese, Spanish and 
French colonies in Africa, the most probable origin of 
T. vivax brought to South America, we have no evidence 
about the route or the livestock species that carried this 
trypanosome to the Americas. The shipment of livestock 
to the New World begin -500 years ago, with inter- 
mediate ports in the Atlantic Islands of Cape Verde and 
Canarias [38-40] allowing the mixture of animals from 
Africa and Europe, where there are no past or present 
day reports of T, vivax. Therefore, the exact geographic 
origin in West Africa, how many times and locations 
of introduction, and the livestock species that carried 
T, vivax to the New World all require further broader 
analysis across Africa and the Americas. 

The population structure observed in the present 
study supports the hypothesis of clonal propagation of 
T, vivax in South America. Evidence is provided by sig- 
nificant levels of linkage disequilibrium between most 
MST loci, the absence of recombinant genotypes and an 
excess of heterozygosity. A clonal structure was also sug- 
gested by a previous MLG analysis from The Gambia, 
where T, vivax was shown to be clonally propagated 
among donkeys [50]. It would be interesting to evaluate 
population structures from other regions of Africa, par- 
ticularly in East African localities where tsetse flies are 
abundant and greater levels of T, vivax diversity occur 
[17,19,34,41,42,59], as well as in endemic tsetse-free areas 
of Sudan and Ethiopia [8-10,59]. 

According to the NJ dendrogram, PCA and STRUCTURE 
analyses, the T, vivax isolates within the main clusters 
were more similar in geographical origin than in date of 
sampling, species of origin, or clinical manifestations. 
The exception was the two EA (Mozambique) isolates 
that clustered together with those from WA. These iso- 
lates also clustered with those from WA and SA using the 
conserved gGAPDH sequences [41,42]. The phylogenetic 
analysis based on Proline-Racemase (TviPRAC) gene se- 
quences showed the two EA isolates included in this study 
clustering closer to SA and WA genotypes than to other 
EA isolates (Caballero et al in preparation). Further 
evidence that some EA isolates are closely related to 
WA/SA genotypes came from a recent study of T, vivax in 
Ethiopia showing isolates from tsetse-free areas sharing 
conserved TviPRAC sequences with WA/SA isolates, 
whereas isolates from tsetse-infested regions showed diver- 
gent sequences [59]. Therefore, increasing evidences are 
supporting the existence in EA of T, vivax genotypes ran- 
ging from very closely related to highly distant to WA/SA 
populations. In addition, results based on polymorphic ITS 
rDNA sequences were sufficiently polymorphic, in agree- 
ment with unique allelic composition, to clearly distin- 
guish the two EA samples included in this study [17,19,34]. 



There is no evidence that the highly divergent T. vivax 
that has been described from East Africa was introduced 
into the South America. The greater genetic diversity of 
T, vivax in East Africa compared to West Africa may be 
related to natural transmission cycles involving a range 
of wild ungulate, and tsetse species in natural reserves of 
wildlife [17,19,35,41,42]. Further analyses are required to 
evaluate whether comparable diversity also occurs in EA 
livestock and to evaluate whether diversity in West Africa 
has been underestimated because all isolates examined 
were from livestock. Studies across African countries are 
of fundamental importance to better understand the 
events shaping the genetic repertoire of T, vivax over a 
long evolutionary history in Africa Similarly, more com- 
prehensive analyses are required to explore the genetic 
diversity of T, vivax across South and Central America to 
better hypothesize about the history of T, vivax in the 
New World. 

Brazilian and Venezuelan populations shared a much 
more similar allelic composition compared to those found 
in West Africa, providing evidence that the T, vivax popu- 
lation we have investigated in South America originated 
from WA genotypes of a common ancestry. Recent com- 
mon ancestors is also suggested by the fact that all SA 
MLGs we have identified can be explained by muta- 
tions that, in general, produced only two repeat MST 
units. This hypothesis is supported by the close rela- 
tionships among T, vivax from the Llanos of eastern 
Venezuela to a vast geographical range including north- 
ern, central, north-eastern, south-eastern, and southern 
Brazilian regions. Given the significant genetic diversity 
found in WA compared to SA countries, the introduction 
of divergent genotypes would most likely have resulted in 
greater genetic diversity in SA T, vivax populations, unless 
there was strong selection for particular genotypes post- 
importation. The more homogeneous SA population 
of T, vivax could result from a population bottleneck 
effect when the parasite adapted to exclusively mech- 
anical transmission. 

The comparison of T, vivax from asymptomatically to 
fatally infected livestock from endemic settings and out- 
breaks, respectively, revealed that of the 6 MLGs de- 
tected in Brazil, only three (MLGs 1-3 sharing highly 
similar allelic composition) were responsible for the out- 
breaks. This result could suggest that outbreaks can occur 
due to clonal expansions underlying more virulent popu- 
lations. However, these three MLGs were also found in 
asymptomatic animals from endemic areas. The sharing 
of MLG between endemic and outbreak sites supports the 
hypothesis that the parasite sources of the outbreaks were 
asymptomatic animals from endemic settings introduced 
into r. vivax-free regions causing severe acute disease, re- 
gardless of the genotypes introduced, into the naive hosts. 
Although our study found no strong association between 
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disease and genotype, the analyses of additional loci and 
samples may reveal genotypes varying in virulence and 
pathogenicity. In fact, sheep and calves experimentally 
infected with a MLGl Brazilian isolate from the Pantanal 
region exhibited very low parasitemia and lacked signs 
of pathogenicity [31]. In contrast, MLG2 isolates from 
cattle and sheep from outbreaks in the semiarid re- 
gion have been shown to be highly virulent for sheep 
and goats [12,13,27,28]. 

It will be useful to assess the temporal stability of T. 
vivax genotypes in longitudinal surveys in specific hosts 
and geographic areas. Our preliminary data did not indi- 
cate substantial temporal changes of genetic profiles of 
parasites in livestock living in areas of enzootic stability. 
Indeed, the highly prevalent MLGl genotype detected in 
Amazonian buffaloes in 2008 and 2009 was also found 
ten years earlier in cattle from the Pantanal, which is 
also an endemic area, and in cattle from an outbreak of 
acute disease with nervous compromises in the south- 
eastern region that occurred in 2008 [29,31]. In contrast, 
in a farm in non-endemic north-eastern Brazil, where 
MLG2 was isolated from an outbreak of high mortality in 
sheep in 2008, only MLGs 4 and 5 (genotypes never found 
in endemic areas) were found in asymptomatic sheep, buf- 
falo and cattle one year after the outbreak, when only the 
sick animals were submitted to treatment [20]. 

Despite our efforts to obtain a representative number 
of isolates, limitations of this study were the very low 
parasitemia and the small number of isolates and loci 
examined. In addition, besides the low number of WA 
samples, they were limited to cattle and results could be 
biased towards isolates selected by these animals. Host 
selection has been considered as an important determin- 
ant of the population structure of T, brucei [60]. Never- 
theless, this is the first MLG analysis comparing SA 
populations, and the results are the first step towards 
the understanding of population structure and genotype 
repertoire of T, vivax throughout distinct epidemiological 
scenarios. However, more substantial sampling from spe- 
cific host species and regions are needed to avoid compli- 
cations of sub-structuring by host, time or space. The 
recent definition of unique repertoires of VSG [61], trans- 
sialidase enzymes involved in parasite virulence [25], and 
genes differentially expressed between T, vivax from 
Venezuela and T, vivax Y486 from West Africa [62] have 
all provided new opportunities to select new markers that 
are useful for comparing the phenotypic and genotypic di- 
versity of T, vivax. 

Conclusions 

Knowledge of the genetic structure of populations is 
critical to investigate the origin, dispersion and impact 
of genetic variation on pathogenicity and epidemiology 
of T, vivax in South America. The MLGs analysis from 



T. vivax across Brazil and Venezuela support both clonal 
propagation and the hypothesis that the isolates exam- 
ined here derived from highly closely related ancestors 
recently introduced from West Africa into the Americas. 
Genetic repertoire was reduced in South America com- 
pared to West Africa. Here, we conducted the first mo- 
lecular comparison of T. vivax from asymptomatically 
to fatally infected livestock that exhibited a range of 
haematological and neurological disorders. Our findings 
found no strong association between genotype, host spe- 
cies, virulence and pathogenicity. 

The increasing number of Brazilian T. vivax outbreaks 
resulting in high mortality of cows, sheep, goats and 
horses, highlights the importance of adopting approaches 
to monitor the spread of T. vivax and the possible selec- 
tion and emergence of genotypes. The molecular markers 
employed in this study are valuable for assessing the 
genetic diversity of American and African populations, 
for reconstructing the pathways of T. vivax introduction 
and dispersion into the Americas, and for determin- 
ing whether particular genotypes emerged locally or 
were imported allowing tracking the source of para- 
sites in outbreaks. 
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